AO-AlSo  187 
UNCLASSIFIED 


MITRE  CORP  MCLEAN  V A  METREK  DIV  F/e  17/7 

CONFLICT  MONITORING  ANALYSIS  OF  PARALLEL  OPPOSITE  DIRECTION  ROU--ETC(U) 
AUG  82  A  P  SMITH  •  DTFA01-82-C-10003 

MTR-82WOOH4-VOL-2  FAA-EM-B2-23-V0L-2  N1 


Afort  No.  FAA-flH-82-23 

VrinsK 


CONFUCT  MONITORING  ANALYSIS 

OF 

PARALLEL  OPPOSITE  DIRECTION  ROUTES 


Arthur  P.  Smith,  III 


The  MURE  Corporation 
McLean,  Virginia  22102 


AUGUST  1982 


Document  is  available  to  the  U .  S .  public  through 
the  National  Technical  Information  Service 
Springfield,  Virginia  22161 


Prepared  for 

U&  DEPARTMENT  OF  TRANSPORTATION 
HDERAL  AVIATION 

OTOE  OF  SYSTMS  EN6MBMN6  MANAGEMENT 
Wastagtan,  DJC.  20591 


*>2  I 


OCT  121982 


1  1  Q 


d.  G«vtrMIMf  ACCMIIM  Ni. 


P  AtJ'O  ft  1 


4.  Ti  tla  an4  Swb>i«l« 


Conflict  Monitoring  Analysis  of  Parallel  Opposite 
Direction  Routes 


/.  Avftitr's) 

Arthur  P.  Smith,  III 


9.  P#cf*npin§  Orf«nit«tiwi  N«n«  m4  A Mr«»s 

The  MITRE  Corporation 
Me trek  Division 
1820  Dolley  Madison  Boulevard 
McLean,  Virginia  22102 


12.  Spwxwint  Afancy  Nam  and  AMrtu 

Office  of  Systems  Engineering  Management 
Federal  Aviation  Administration 
Department  of  Transportation 


Tcciuiicsl  Repeat  Deceewetetion  Pope 


3.  flacryiant't  Catalaf  Na. 


•rf«nRin|  OrfcniMfiwi  C»4« 


9.  OrfifiixatiM  N». 

MTR-82W00114 ,  Volume  II 


10.  Unit  N#.  (TRAI5) 


1 1 .  Guttract  or  Grant  N«. 

DTFA01-82-C10003 


*9.  T yp0  of  Rooort  mnd  Poriod  C ovorod 


IS.  SwyyUwnfry  N«t«a 


1 6.  Abstract 

This  paper  reports  on  the  development  of  the  Conflict  Monitoring  Analysis.  A  previous 
report  (FAA-EM-80-16)  described  the  estimates  of  the  probability  of  horizontal  overlap 
and  controller  intervention  rate  for  same  direction  adjacent  parallel  routes.  This  ! 

work  extends  that  methodology  to  opposite  direction  adjacent  parallel  routes.  For  j 

both  the  probability  of  horizontal  overlap  and  the  controller  intervention  rate, 
trial  results  based  on  data  are  given. 
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1.  INTRODUCTION 


The  first  volume  of  this  report  presents  results  of  the  work  to 
extend  the  Conflict  Monitoring  Analysis  to  opposite  direction 
traffic  flows  at  the  same  flight  levels  for  parallel  routes. 

The  appendices  in  this  volume  present  the  details  of  the 
opposite  direction  analysis.  As  pointed  out  in  Volume  I,  much 
of  the  supporting  analysis  for  the  opposite  direction  study  is 
common  to  the  earlier  study  of  same  direction  routes  and  can  be 
found  in  Volume  II  of  Reference  3.  The  appendices  in  this 
volume  address  the  analytical  differences  between  same  and 
opposite  direction  traffic  flows.  The  specific  topics  discussed 
in  this  volume  are  described  below: 

o  Some  initial  work  to  address  conflicts  involving 

turning  aircraft  was  done  in  the  opposite  direction 
analysis.  Appendix  A  describes  the  accelerations 
observed  in  the  data. 

o  In  the  same  direction  analysis,  errors  in  the 

estimates  of  the  crcnstrack  separations  and  closing 
speeds  were  considered.  In  the  opposite  direction 
analysis,  the  alongtrack  dimension  also  becomes 
important  since  the  closing  speed  is  largely  in  that 
dimension.  Appendix  B  addresses  the  surveillance 
system  and  tracker  errors  which  contribute  to  the 
errors  in  the  estimates  of  the  crosstrack  and 
alongtrack  separation  and  closing  speeds. 

o  By  considering  both  the  crosstrack  and  alongtrack 
dimensions  in  the  opposite  direction  analysis,  the 
description  of  the  conflict  region  is  also  modified 
from  the  same  direction  analysis.  Appendix  C 
describes  the  opposite  direction  conflict  surface.  As 
in  the  same  direction  analysis  the  delay  incurred 
before  the  controller  realizes  that  the  aircraft  pair 
is  within  the  conflict  region  is  required.  Appendix  D 
addresses  the  computation  of  the  delay  for  the 
opposite  direction  conflict  region. 

o  In  the  same  direction  analysis  the  aircraft  performing 
the  resolution  maneuver  was  assumed  to  turn  toward  its 
assigned  route  centerline.  In  the  high  closure  rate 
situation  of  opposite  direction  traffic  flow  this 
maneuver  may  not  be  the  most  effective.  Appendix  E 
desrlbes  the  Conflict  Resolution  Advisory  which  is 
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being  developed  for  NAS  and  Is  used  In  the  opposite 
direction  Conflict  Monitoring  Analysis  as  the  logic 
for  selecting  the  resolution  maneuver. 

o  The  net  result  of  the  Conflict  Monitoring  Analysis  is 
to  estimate  the  probability  of  horizontal  overlap. 

The  inclusion  of  acceleration  and  both  left  and  right 
turns  for  resolutions  means  that  there  are  new 
descriptions  of  the  overlap  regions  in  the  opposite 
direction  analysis.  These  are  discussed  in  Appendix  F. 

o  The  combination  of  all  the  appropriate  probabilities 
to  arrive  at  the  estimate  of  the  probability  of 
horizontal  overlap  is  addressed  in  Appendix  G. 

An  indication  of  the  parts  of  the  opposite  direction  analysis 
that  use  the  results  from  these  appendices  is  given  in  Figure 
1-1. 
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APPENDIX  A 


AVERAGE  RELATIVE  LATERAL  ACCELERATIONS 


Even  though  the  opposite  direction  Conflict  Monitoring  Analysis  is 
not  developed  to  the  point  of  making  satisf icatory  estimates  of 
horizontal  overlap  probability  with  accelerated  flight  trajectories, 
this  appendix  will  document  the  accelerations  found  in  one  sample  of 
flight  tracks.  In  particular  it  is  of  Interest  to  examine  the 
relative  lateral  (crosstrack)  accelerations  averaged  over  a  two 
minute  period.  The  averaging  is  desired  because  the  conflict 
encounter  lasts  about  two  minutes  or  until  the  aircraft  pass. 
Furthermore,  an  average  acceleration  estimate  is  a  more  stable 
estimate  than  an  instantaneous  acceleration  estimate  made  from 
position  reports. 

The  average  acceleration  estimates  were  made  by  considering  184 
aircraft  tracks  on  jet  route  J146  in  Cleveland.  The  data  available 
were  the  crosstrack  deviations  from  the  route  centerline  at  twelve 
second  time  intervals.  This  set  of  points  was  smoothed  with  a  cubic 
spline  function.  Between  data  points  1  and  i+1  the  spline  function 
is  defined  as 

yi  -  Ai  +  B^(t-tj)  +  Cj^t-ti)2  +  Dj/t-tj)-* 

tilt<t1+1  (A-l) 

such  that  y  and  its  first  and  second  derivatives  are  continuous  at 

each  point.  The  technique  used  here  is  called  spline  function 
smoothing  (Reference  2)  where  the  objective  is  to  minimize  the 
second  derivative  of  the  function  over  the  entire  interval. 

Given  the  crosstrack  deviation  between  each  pair  of  points  given  by 
(A-l)  the  average  crosstrack  acceleration  between  point  i  and  point 
i+1  is 

®i  “  2ci  +  3DI  (ci+l  “  ci)  (A-2) 

In  this  case  tj+i  -  tj-  is  equal  to  12  seconds.  After  the 

average  lateral  acceleration  is  computed  over  each  12  second 
interval  an  average  over  ten  12  second  intervals  (2  minutes)  is 
taken.  After  a  gap  of  two  minutes  in  the  track  another  two  minute 
average  is  taken.  The  resulting  histogram  of  average  lateral 
accelerations  for  single  aircraft  is  shown  in  Figure  A-l.  Although 
the  samples  in  the  histogram  are  not  truly  independent  they  will  be 
considered  as  such  for  the  purposes  of  this  analysis. 
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Since  Che  relative  lateral  acceleration  is  required,  the  histogram 
in  Figure  A-l  was  convolved  with  itself  using  the  Fast  Fourier 
Transform  (see  Appendix  0  of  Volume  II  of  Reference  3).  The 
resulting  histogram  which  was  presented  in  Volume  I  of  this  report 
is  shown  in  Figure  A-2. 


AVERAGE  RELATIVE  CROSSTRACK 
ACCELERATION  HISTOGRAM 


APPENDIX  B 


SURVEILLANCE/TRACKER  ERRORS 

B.l  Background 

In  the  same  direction  Conflict  Monitoring  Analysis  we  needed  to 
know  the  probability  of  observing  the  aircraft  pair  outside  the 
conflict  boundary  when  the  aircraft  pair  was  truly  inside  the 
boundary.  The  uncertainty  of  observing  the  aircraft  pair  inside 
or  outside  the  conflict  boundary  is  due  to  the  errors  in  the 
surveillance/tracking  system.  It  was  assumed  in  the  same 
direction  case  that  the  errors  in  the  crosstrack  closing  speed 
and  the  crosstrack  separation  were  sufficient  to  describe  the 
situation  .  Therefore  the  estimate  of  the  probability  of 
observing  the  aircraft  pair  outside  the  conflict  boundary  when 
it  was  truly  inside  was  made  using  an  approximation  for  the 
volume  under  a  correlated  bivariate  normal  distribution.  The 
procedure  for  making  that  estimate  is  discussed  in  Appendix  C  of 
Volume  II  of  Reference  3. 

In  the  opposite  direction  Conflict  Monitoring  Analysis  there  is 
an  analogous  need  for  Che  probability  of  being  observed  outside 
the  conflict  boundary  when  the  aircraft  pair  is  truly  inside  the 
boundary.  However,  in  the  opposite  direction  case  the  errors  in 
the  alongtrack  separation  and  the  crosstrack  speed  of  one  of  the 
aircraft  due  to  the  surveillance/tracking  system  are  also 
important.  Therefore  the  probability  of  observing  a  conflicting 
pair  outside  the  conflict  surface  will  be  estimated  using  a 
correlated  quadrivariate  normal  distribution. 

This  appendix  will  discuss  how  the  parameters  of  the 
quadrivariate  normal  distribution  were  arrived  at.  This  is 
followed  by  a  discussion  of  how  one  can  randomly  sample  from  a 
correlated  quadrivariate  normal  distribution. 

B.2  Parameter  Estimation 


The  parameter  of  the  correlated  quadrivariate  normal 
distribution  were  estimated  via  simulation.  The  simulation  was 
run  for  four  orientations  of  the  parallel  routes  with  respect  to 
the  radar.  In  each  case  15  pairs  of  aircraft  were  flown  down 
the  routes.  The  aircraft  tracks  came  from  the  FAA's  VOR 
navigation  data  base.  There  were  200  randomly  chosen  tracks 
from  the  selected  routes  in  the  Cleveland  ARTCC.  From  two 
opposite  direction  routes  50  aircraft  were  chosen  to  comprise 


Che  population  (approximately  25  airc  sft  on  each  route).  In 
the  simulation  the  pairs  of  aircraft  were  chosen  from  their 
respective  route  populations.  These  pairs  were  chosen  with 
replacement.  The  aircraft  were  placed  on  their  respective 
routes  within  5  nmi  alongtrack  of  the  designated  starting  points 
of  the  routes.  The  aircraft  were  allowed  to  travel 
approximately  10  nmi  before  statistics  were  taken.  This  allowed 
for  the  transients  in  the  tracker  to  die  out.  Figure  B-l  shows 
the  configurations  of  the  four  cases. 

The  errors  in  the  surveillance  system  were  modelled  in  the 
following  manner.  It  was  assumed  that  the  azimuth  error  was 
normally  distributed,  mean  of  zero  and  standard  deviation  of 
0.264  degrees  (3  ACPs).  The  range  error  is  a  combination  of  two 
errors.  The  first  is  a  bias  error  built  into  the  airborne 
transponder.  This  is  assumed  to  be  uniformly  distributed 
between  ♦  0.5  of  a  microsecond  (  +  0.08  nmi).  This  error  is 
assigned  to  each  aircraft  as  it  enters  the  simulation.  The 
second  range  error  is  due  to  the  range  quantization  in  the 
common  digitizer.  The  quantization  is  1/8  nmi  and  the  error  is 
applied  at  each  radar  update.  The  probability  of  a  radar  return 
(blip/scan)  is  set  at  0.95. 

Each  track  contains  the  information  about  the  position  and 
velocity  of  the  aircraft  at  each  update  time.  At  each  update 
the  "true"  position  and  velocity  components  are  compared  to  the 
"observed"  position  and  velocity  components.  Positive  values 
for  the  variables  have  the  following  meanings: 

Alongtrack  separation  (x)  —  the  aircraft  appear  to  be 
farther  apart  than  they  truly  are, 

Croastrack  separation  (y)  —  the  aircraft  appear  to  be 
farther  apart  than  they  truly  are, 

Croastrack  closing  speed  (y)  —  the  aircraft  appear  to  be 
closing  faster  than  they  truly  are, 

Croastrack  speed  of  aircraft  1  (y^)  —  the  aircraft 
appears  to  be  flying  north  faster  than  it  truly  is. 

The  results  for  the  sixty  pairs  of  aircraft  (15  pairs  times  4 
route  orientations)  are  given  in  Tables  B-l,  8-2,  and  B-3. 

UX  is  the  estimate  of  the  mean  alongtrack  separation  error, 

UY  is  the  estimate  of  the  mean  crosstrack  separation  error, 
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UYD  is  the  estimate  of  the  mean  crosstrack  closing  speed 
error, 

UYD1  is  the  estimate  of  the  mean  crosstrack  speed  error  of 
aircraft  one, 

SX,  SY,  SYD,  SYD1  are  the  estimates  of  the  corresponding 
one  sigma  values, 

PXY  is  the  estimate  of  the  correlation  coefficient  between 
the  alongtrack  and  crosstrack  separation  errors, 

PXYD  is  the  estimate  of  the  correlation  coefficient  between 
the  alongtrack  separation  and  the  crosatrack  closing  speed 
errors, 

PXYD1  is  the  estimate  of  the  correlation  coefficient 
between  the  alongtrack  separation  and  the  crosstrack  speed 
of  aircraft  one  errors, 

PYYD  is  the  estimate  of  the  correlation  coefficient  between 
the  crosstrack  separation  and  the  crosstrack  closing  speed 
errors, 

PYYD1  is  the  estimate  of  the  correlation  coefficient 
between  the  crosstrack  separation  and  the  crosstrack  speed 
of  aircraft  one  errors, 

PYDYD1  is  the  estimate  of  the  correlation  coefficient 
between  the  crosstrack  closing  speed  and  the  crosstrack 
speed  of  aircraft  one  errors. 

Table  B-l  shows  the  mean  and  standard  deviations  of  the  mean 
errors.  Assuming  the  means  are  distributed  normally,  a  t-test 
could  be  used  to  test  the  hypothesis  that  the  means  have  been 
sampled  from  a  population  with  a  zero  mean.  Even  though  two  of 
the  t-tests  would  statistically  reject  the  hypothesis  of  zero 
means,  operationally  the  means  are  near  enough  to  zero  that  a 
value  of  zero  will  be  used  in  the  analysis. 

Table  B-2  shows  the  means  and  standard  deviations  of  the 
standard  deviation  of  the  errors.  For  both  the  alongtrack  and 
crosstrack  separation  errors  the  standard  deviation  is  estimated 
to  be  about  0,3  nml.  The  standard  deviation  of  the  crosstrack 
closing  speed  error  is  about  39  kts«  The  standard  deviation  of 
the  crosstrack  speed  of  aircraft  one  is  about  23  kts. 
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Table  B-3  ahovs  Che  means  and  standard  deviations  of  the 
pairwise  correlation  coefficients.  If  normality  is  assumed,  a 
t-test  could  be  used  to  test  the  hypothesis  that  the  crosstrack 
errors  are  uncorrelated.  This  test  shows  that  the  alongtrack 
dimension  is  basically  uncorrelated  with  the  crosstrack 
dimension.  This  is  reasonable  since  the  computation  of  the 
crosstrack  velocities  is  based  only  on  the  crosstrack  positions 
and  the  crosstrack  errors  are  uncorrelated  with  the  alongtrack 
errors. 

Thus,  for  the  purposes  of  the  Conflict  Monitoring  Analysis,  the 
quadrivarlate  normal  distribution  parameters  listed  in  Table  B-4 
are  used  to  specify  the  surveillance/tracker  errors  for  a  pair 
of  opposite  direction  aircraft. 

B.3  Sampling  From  a  Quadrivarlate  Normal  Distribution 

In  the  determination  of  detection  delays  (Appendix  D  of  this 
paper)  it  will  be  necessary  to  sample  from  a  correlated 
quadrivarlate  normal  distribution.  The  reason  this  is  necessary 
is  that  there  is  no  general  approximation  for  a  volume  under 
such  a  distribution. 

If  U^,  U2  Um  are  Independent  unit  normal  variables 

their  joint  density  is 


eu(u)-(2K)-a/2>“  «,[-  5  £  0  ]  ] 


“(2  7r  )-(1/2)  exp  (  “  |  UTJi  )  . 


If  we  apply  a  nonsingular  linear  transformation^  to 
Yt  -  <Ylt...,  Ym), 

UT  -  YTHT  with  I H |  t  0,  (B-l) 
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TABLE  B-4 


PARAHETER  VALDES 
OP  THE 

QOADRIVABIATE  SOB VEILLAHCE/TH ACKER 
ERROR  DISTRIBUTION 


f — 

r  — 

1 

* 

Pa rase ter 
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1- 

1 

J 

Kx 

1 

1 

I  0  ui 

1 
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i 
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1 

. 1 

r 

a* 

1 

|  0.3  mi 

1 

<ry 

|  0.  3  mi 

at 

|  39  kts 

'*1 

I  23  kts 

1  -  • 

1 

|  _ 

^  -r| 

r“ 

1 

1  0 

*1 
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1  o 

1  o 

PT7 

t  -0.6 

Pyf, 

I  -0.4 

1  / Route  1\ 

|  (North  of  j 

I  \  Route  2/ 

p*  • 

YYl 

1  0.5 

I  /  Route  1\ 

1  (forth  of ) 

I  \  Route  2/ 

i _ 

1 

-X. 

_ i 

B-9 


we  find  that  Y  has  a  joint  density  function 


P (y)  “  (2tt)  exp(-7yTHTH^) 

-(-)m  -  T 

=  (2n)  | R |  'exp(-7 yT R y ) 

where  R**HTH  is  positive  definite. 

The  variance-covariance  matrix  of  Y  is  R”*.  Since  the 
correlation  coefficient  is  related  directly  to  the  covariance 

.  SUEZ _ 

Pxy  <rx  try 

We  see  that  a  matrix  of  the  form 


1 


P 


Pml 


can  represent  the  symmetry  of  R-1.  Thus,  knowing  P  we  can 
invert  to  find  R.  Knowing  R  we  solve  g  *  H^H  to  find  H. 

Then  knowing  g  we  will  know  the  transformation  to  convert  the 
independent  unit  normal  variables  (U)  to  the  correlated 
variables  (Y): 


UT(HT)'1-YT. 


Inverting  P  to  find  R  is  straightforward.  The  solution  of  (B-2) 
is  accomplished  by  a  Cholesky  decomposition  which  is  a  matrix 
manipulation  supported  by  "SAS"  (the  Statistical  Analysis  System 
computer  program  package) . 

As  an  example  of  the  generation  of  a  correlated  quadrivariate 
random  numbers  consider  four  variables:  x,  y,  y,  and  y^. 

Assume  that  x  is  independent  of  the  other  variables. 

Let  the  P  matrix  be 


‘1 

0 

0 

0  -| 

0 

1 

-0.66 

0.59 

0 

-0.66 

1 

-0.67 

.0 

0.59 

-0.67 

1 
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The  resulting  matrix  (HT)_1  is 


’I 

0 

0 

0" 

0 

0.724400 

0 

0 

(hV1  - 

0 

-0.356565 

0.742361 

0 

This  means  that 

.0 

0.590000 

-0.670000 

1. 

Ul*x 

(0.7244  U2  -  0.356565  U3  +  0.59  U4)  try 
(0.742361U3  -  0.67  U4  )  0^ 


[hm  »4ffh 


(B-3) 


where  U.  are  unit  normal  random  numbers  and  crj  are  the 
standard  deviations  of  each  variable.  If  crx  “1,  cry  “  2, 
a-  ■  3,  and  -  4  and  a  Monte  Carlo  simulation  is  run  150 

tlsaes  picking  x, 
are:  1 


Variable 

x 

y 

y 

*1 


y,  and  ylf  according  to  (B-3)  the  result  1 

Mean 

Standard  Deviation 

-0.0009  (0) 

0.98  (1) 

0.0741  (0) 

2.22  (2) 

-0.4276  (0) 

3.01  (3) 

0.5819  (0) 

4.00  (4) 

Correlation  Coefficents 


xy 

—0.022 

(0) 

xy 

—0.033 

(0) 

**1 

-  0.087 

(0) 

yy 

—0.709 

(-0 

yfi 

-  0.609 

(  0 

yyi 

—0.678 

(-0 

when  the  values  in  parentheses  are  the  parameters  of  the 
distribution  being  sampled  from. 
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APPENDIX  C 


THE  CONFLICT  SURFACE 


This  appendix  will  develop  the  analysis  necessary  to  describe  the 
conflict  surface  for  opposite  direction  traffic.  The  conflict 
surface  is  defined  by  those  values  of  the  separation  and  closing 
speed  between  two  aircraft  which,  when  rectilinearly  projected  ahead 
in  time,  would  Indicate  that  the  aircraft  would  be  separated  by  5 
noi  in  2  minutes  or  less.  Just  outside  such  a  conflict  surface  the 
aircraft  would  not  be  projected  to  be  separated  by  S  nml  within  2 
minutes.  To  understand  the  conflict  surface  it  is  necessary  to 
define  the  sign  conventions  of  the  parameters  and  explictly  describe 
the  conflict  criteria. 


C.l  Sign  Conventions 


The  following  sign  conventions  are  used  in  the  analysis: 

y^  crosstrack  deviation  of  aircraft  1.  Positive  is 

away  from  own  route  toward  the  other  route. 

y2  crosstrack  deviation  of  aircraft  2.  Positive  is 

away  from  own  route  toward  the  other  route. 

y  crosstrack  separation  ■  RS-Cy^  +  y2)  where  RS 

is  the  route  spacing.  Positive  indicates  that 
the  aircraft  are  oriented  with  respect  to  each 
other  in  the  same  way  their  respective  routes  are 
oriented.  In  other  words,  if  route  2  is  north  of 
route  1  and  aircraft  2  is  north  of  aircraft  1, 
then  y  is  positive. 

yj  is  the  crosstrack  velocity  of  aircraft  1. 

Positive  is  away  from  own  route  toward  the  other 
route . 

y2  is  the  crosstrack  velocity  of  aircraft  2. 

Positive  is  away  from  own  route  toward  the  other 
route . 


y-yi+72 


is  the  "closing  speed”  of  the  aircraft,  y  is 
positive  if  aircraft  are  closing  and  y  >  0. 
y  is  positive  if  aircraft  are  opening  and  y  <  0 


X}  is  the  alongtrack  position  of  aircraft  1. 

Positive  is  in  direction  of  flight  of  aircraft 
1.  The  origin  is  arbitrary. 

X2  is  the  alongtrack  position  of  aircraft  2.  The 

reference  frame  is  the  same  as  for  aircraft  1. 

x*X2~xl  Is  the  alongtrack  separation.  Positive  indicates 

aircraft  have  not  yet  passed. 

x^  is  alongtrack  velocity  of  aircraft  1  (+) 

X2  is  alongtrack  velocity  of  aircraft  2  (-  if 

aircraft  2  is  flying  in  a  direction  opposite  that 
of  aircraft  1.) 

x«*X2-xl  is  the  alongtrack  "closing"  velocity.  Minus 

means  closing  if  X2>xj. 

C.2  Conflict  Criteria 

A  pair  of  aircraft  is  said  to  be  in  potential  conflict  if  the 
two  aircraft  are  projected  to  be  within  5  nmi  (D)  within  2 
minutes  (T).  If  one  projects  the  aircraft  assuming  rectilinear 
motion,  one  can  plot  the  aircraft  separation  as  a  function  of 
time  as  shown  in  Figure  C-l.  If  the  projected  trajectory  of  the 
pair  goes  through  the  shaded  area  the  pair  would  be  in  potential 
conflict  (i.e.,  at  some  point  in  time  less  than  T  the  separation 
of  the  aircraft  would  be  less  than  D).  There  is  an  obvious 
distinction  that  one  could  make  between  various  trajectories  — 
those  that  do  and  those  that  don't  pass  through  the  shaded 
area.  There  are  three  types  of  trajectories  that  do  pass 
through  the  shaded  area.  These  are  characterized  by  the  time  of 
closest  approach  (minimum  d)  of  the  trajectory.  As  shown  in 
Figure  C-l,  the  time  of  closest  approach  for  trajectory  1  is 
greater  than  T.  Trajectory  2's  time  to  closest  approach  is  less 
than  T  and  the  separation  of  the  aircraft  at  T  is  less  than  D. 
Trajectory  3's  time  to  closest  approach  is  less  than  T  but  the 
separation  la  greater  than  D  at  time  T.  In  any  case,  if  there 
is  a  potential  conflict  the  distance  of  closest  approach  is  less 
than  0.  The  converse,  however,  is  not  true.  A  distance  of 
closest  approach  less  than  D  does  not  necessarily  mean  a 
potential  conflict. 


r 


There  are  three  events  that  one  has  to  consider  when  determining 
whether  a  trajectory  is  in  potential  conflict.  One  event  is 
"the  time  to  closest  approach  is  less  than  T,"  (tca<T).  The 
second  event  is  "the  distance  of  closest  approach  is  less  than 
D,”  (dca<D).  The  third  event  is  “the  distance  at  time  T  is 
less  than  D,”  (dT<D).  The  conflict  condition  is  thus: 

[(tc.i<  T)  n  (dca<D)]u  (dT<  D).  (C-l) 

This  logical  expression  will  be  used  to  determine  the  conflict 
region  boundary  in  the  following  analysis. 

C.3  Development  of  the  Conflict  Boundary 

In  the  development  of  the  conflict  boundary  it  is  assumed  that 
the  forward  speeds  of  each  aircraft  are  constant  and  not 
necessarily  equal.  Thus,  by  specifying  Sj.,  S2,  (the  two 
forward  speeds),  y,  yj_,  x,  and  y  the  complete  geometry  of  the 
encounter  is  known.  This  means  that  x  is  a  function  of  Sj, 

S  ,  y,  and  y^  .  Thus,  given  S^,  S2,  x,  y^  the 
conflict  geometry  is  described  by  y  and  y. 

If  the  aircraft  are  separated  at  t  »  0  by  x,  y  then  the  aircraft 
will  be  separated  by 

d  -  7 (x  +  xt)2  +  (y  -  yt)2  (C-2) 

at  time  t.  This  assumes  rectilinear  motion.  From  (C-2)  one  can 
compute  the  time  to  closest  approach.  It  is 

cca  •  (yy  -  x*)  /  (*2  +  y2)  (c-3) 

If  we  let  tca  *  T  and  solve  (C-3)  for  y  we  get 

y  -  (T(x2  +  y2)  +  xx)  /  y  (C-4) 

Equation  (C-4)  specifies  the  locus  of  points  where  the  time  to 
closest  approach  is  exactly  T  minutes  ahead.  Figure  C-2  shows 
the  two  configurations  of  this  locus  of  points  for  the  given 
values  of  S^,  S2,  x,  and  y^.  The  shaded  areas  represent 
those  combinations  of  y  and  y  where  tca<  T.  Being  in  a  shaded 
region  in  Figure  C-2  is  one  of  the  events  in  (C-l)  that 
determine  whether  the  aircraft  pair  will  be  in  potential 
conflict . 


The  second  condition  that  one  has  to  consider  when  determining 
if  a  potential  conflict  exists  is  if  the  closest  approach 


400  kts 
600  kts 
100  kts 


distance  is  less  than  D.  This  condition  is  mathematically 
expressed  by  substituting  tca  as  given  by  (C— 3)  for  t  in  (C-2) 
and  letting  d  in  (C-2)  equal  D.  After  rearranging  terms  we 
arrive  at 

y  -  (-xy  +  Dv/x^  +  y2)/x  (C-5) 

Thus,  equation  (C-5)  describes  the  loci  of  points  where  the 
distance  of  closest  approach  is  equal  to  D.  Between  the  two 
lines  described  by  (C-5)  in  the  y,  y  plane,  the  closest  approach 
distance  is  less  than  D.  In  terms  used  in  Volume  I  of  this 
report,  this  condition  defines  the  ‘’backside1’  region  of  the 
conflict  surface. 

The  third  condition  specified  in  (C-l)  for  a  potential  conflict 
is  that  the  separation  of  the  aircraft  at  the  look-ahead  time  T 
is  less  than  D.  This  condition  is  mathematically  described  by 
letting  t  ■  T  and  d  ■  D  in  equation  (C-2).  After  solving  for  y 
we  get 

y  -  yT  +  xT)2  (C-6) 

In  the  y,  y  plane  this  expression  describes  two  lines.  The  area 
between  the  two  lines  is  where  the  separation  is  less  than  D  at 
time  T.  This  condition  defines  the  "leading  edge”  region  of  the 
conflict  surface. 

As  an  illustration  of  the  use  of  (C-l),  (C-4),  (C-5),  and  (C-6) 
in  the  determination  of  the  conflict  boundary,  consider  Figure 
C-3.  In  this  figure  the  forward  speeds  of  the  two  aircraft  are 
■  400  kts,  S2  *  600  kts.  Also  x  ■  33  nmi  (the  alongtrack 
separation  )  and  yj^  *  100  kts  (the  crosstrack  speed  of 
aircraft  1).  The  shaded  regions  in  Figure  C-3  defined  by  lines 
4  (corresponding  to  equation  (C-4))  is  where  the  time  to  closest 
approach  is  less  than  T.  The  intersection  of  this  region  with 
the  region  between  lines  5  defines  where  the  distance  of  closest 
approach  is  less  than  D  and  the  time  of  closest  approach  is  less 
than  T.  This  is  shown  by  the  cross  hatched  region  in  Figure  C-3. 

However,  the  complete  conflict  region  is  given  by  the  union  of 
this  region  and  the  region  between  lines  6.  Figure  C-4  shows 
the  complete  conflict  boundary  for  the  given  x  and  yj.  The 
maximum  extent  in  y  for  the  conflict  region  is  determined  by  the 
maximum  value  of  the  single  aircraft  crosstrack  speed,  in  this 
case  200  kts  (200+100  -  300  and  -200+100-100). 
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EXAMPLE  OF  CONFLICT  BOUNDARY 
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It  should  be  obvious  from  equations  (04),  (05),  and  (06)  that 
for  various  values  of  x,  yj_,  y,  and  y  either  equation  (05)  or 
(C-6)  will  determine  the  boundary.  In  general,  as  x  becomes 
small  (less  than  about  28  nml),  lines  6  will  not  exist  because 
the  expression  under  the  radical  in  equation  (06)  is  negative. 
This  means  that  at  exactly  2  minutes  in  the  future  the  aircraft 
will  have  passed  each  other  and  would  be  separated  by  more  than 
5  nml.  In  this  case  the  conflict  boundary  is  determined  by 
equation  (05).  At  larger  values  of  x  (greater  than  about  35 
nml  for  the  given  parameter  values)  there  is  no  intersection 
between  the  regions  where  tca  <  T  and  dca  <  D.  In  this  case 
lines  6  determine  the  conflict  boundary.  At  values  of  x  greater 
than  about  39  nmi  the  conflict  boundary  vanishes  altogether 
because  the  aircraft  are  too  far  apart  to  close  within  D  nmi 
within  T  minutes. 

In  an  attempt  to  graphically  display  the  conflict  surface  Figure 
05  was  produced.  This  figure  shows  the  top  of  the  conflict 
surface  for  one  value  of  y^.  Due  to  the  computer  graphics 
package  used  it  was  not  feasible  to  show  that  the  bottom  of  the 
conflict  surface  Is  symmetrical  with  the  top.  However,  Figure 
05  does  show  the  very  steep  "leading  edge"  of  the  conflict 
surface  (larger  x)  and  the  shallower  "backside"  region  (smaller 
x). 
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APPENDIX  D 


PENETRATION  OF  THE  CONFLICT  SURFACE 

- OTP - 

DETECTION  DELAYS 


This  appendix  will  describe  1)  how  one  can  determine  where  the 
aircraft  pair  trajectory  will  penetrate  the  conflict  surface  and  2) 
how  far  inside  the  conflict  surface  the  pair  can  penetrate  before 
being  detected.  The  basic  assumption  is  that  the  aircraft  pair 
starts  with  a  large  enough  initial  alongtrack  separation  that  it  is 
not  immediately  in  conflict.  Then,  based  on  the  particular  closing 
speeds  and  separations  the  aircraft  pair  trajectory  may  penetrate 
the  conflict  surface.  As  shown  in  Appendix  C  the  conflict  surface 
has  two  distinct  regions.  At  large  alongtrack  separations  there  is 
a  steep  leading  edge  to  the  conflict  surface  in  the  x-y-y  space.  At 
smaller  values  of  alongtrack  (x)  separation  the  conflict  surface  is 
flatter.  This  is  referred  to  as  the  backside  region.  The  aircraft 
pair  trajectory  will  either  1)  miss  the  conflict  surface  altogether, 
2)  penetrate  the  leading  edge  of  the  conflict  surface,  or  3) 
penetrate  the  backside  region.  A  trajectory  without  acceleration 
can  only  penetrate  the  leading  edge  of  the  conflict  surface.  With 
acceleration  the  trajectory  can  penetrate  either  the  leading  edge  or 
the  backside  region.  The  algorithms  for  determining  the  point  of 
penetration  of  the  conflict  surface  are  developed  below. 

D.l  Leading  Edge  Penetration 

This  algorithm  is  defined  using  the  following  parameters:  the 
initial  crosstrack  separation,  the  initial  crosstrack  closing 
speed,  the  crosstrack  closing  acceleration,  the  initial 
alongtrack  separation,  the  initial  alongtrack  closing  speed,  and 
the  alongtrack  deceleration. 

D.1.1  Leading  Edge  Penetration  Without  Lateral  Acceleration 

The  first  step  in  his  algorithm  is  to  assume  no  accelerations  or 
decelerations  and  determine  if  the  trajectory  of  the  aircraft 
pair  penetrates  the  leading  edge  of  the  conflict  surface.  The 
leading  edge  of  the  conflict  surface  is  defined  as 

y  -  yT+  v/  D2  -  (x+xT)2  (D-l) 


where 


y  -  crosstrack  separation 


D-l 


y  ■  crosstrack  closing  speed  (+  is  closing  if  y  >  0) 
x  -  alongtrack  separation 
x  *  alongtrack  closing  speed  (-  is  closing) 

D  -  radar  minimum  separation 
T  »  look  ahead  time 

The  trajectory  without  acceleration  as  a  function  of  time,  t,  is 
giren  as 

y  -  y0  -  y0t  (°“1 2> 

X  -  x0  +  XQt 

where  y0  and  xQ  are  the  initial  crosstrack  and  alongtrack 
closing  speeds  and  x0  and  y0  are  the  Initial  alongtrack  and 
crosstrack  positions,  respectively.  Eliminating  t  from  (D-2)  we 
get 


Substituting  (D-3)  into  (D-l)  and  solving  for  x  we  get 

1  +  <  w* 


(D-3) 


(D-4) 


where 

DISC  -  D2(l+(i/i)2)-((y/i)(j-iT))2 * * 
o  o  o  o  o  o 

-(T-tx  /x  )(2y  -y  T)y  -y  (y  -2y  T) 
v  o  o  Jo  o  7o  O  O  0 

The  conditions  that  must  be  met  for  (D-4)  to  be  a  valid  solution 
for  the  penetration  of  the  leading  edge  are  as  follows: 

1)  DISC  must  be  greater  than  or  equal  to  zero 

2)  The  larger  of  the  two  roots  in  Equation  (D-4)  is  chosen 

3)  The  chosen  root  for  x  must  be  between  zero  and  xQ 

D.1.2  Leading  Edge  Penetration  With  Lateral  Acceleration 


The  value  of  x  found  from  (D-4)  is  used  as  the  starting  value  to 

solve  for  the  penetration  of  the  leading  edge  with 

acceleration.  The  accelerated  trajectory  is 


where  the  variables  are  as  defined  above  with  the  addition  of 

y  -  crosstrack  acceleration  (+  is  closing) 
x  -  a long track  acceleration  (+  is  decelerating) 

The  alongtrack  acceleration  is  approximated  by  assuming  that  the 
crosstrack  track  acceleration  is  partitioned  equally  between  the 
acceleration  will  be  constant.  This  assumption  is  made  to  ease 
the  computational  complexity  and  does  not  introduce  too  great  an 
error  over  the  parameter  ranges  of  interest.  With  these 
assumptions  it  can  be  shown  that 
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(D-6) 


As  the  aircraft  accelerate  the  new  crosstrack  and  alongtrack 
velocities  are  given  by 


x 

y 


x©  +  set 
yQ  + 


(D-7) 


Equations  (D-l),  (D-5),  and  (D-7)  are  used  iteratively  (as 
explained  below)  to  solve  for  t  and  hence  for  x,  y,  x,  and  y. 

If  there  is  no  real  solution  to  (D-4)  or  if  the  accelerations 
are  such  that  no  accelerated  penetration  of  the  leading  edge 
takes  place,  then  the  possibility  of  a  penetration  of  the 
backside  of  the  conflict  surface  is  explored. 

D.2  Backside  Penetration 

The  only  way  for  the  backside  of  the  conflict  surface  to  be 
penetrated  is  with  accelerated  motion.  An  accelerated  motion 
trajectory  is  described  by  equation  (D-5). 


The  backside  of  the  conflict  surface  as  given  in  (C-5)  is 


xy  ■  -xy  +  D  \/x2  +  ^ 
where  x  and  y  are  given  in  (D-7). 


(D~8) 


By  substituting  (D-5)  and  (D-7)  into  (D-8)  and  solving  for  t  one 
gets  the  time  at  which  the  trajectory  penetrates  the  surface. 

The  method  used  to  solve  the  non-linear  equation  Is  that  of 
Newton  where 

cn  ”  cn-l  “  ^^n-l^^'  (*n-l7 
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where 


tn  is  Che  time  at  the  nth  interation 
f  (  )  ■  0  is  the  function  whose  root  is  to  be  found 
f'(  )  is  the  derivative  of  the  function  with  respect  to 
time. 

When  a  value  of  t  is  found  that  satisfies  (D-8)  within  a  stated 
tolerance,  it  is  substituted  into  (D-5)  and  (D-7)  to  give  values 
for  x,  y,  x,  and  y.  Assuming  that  the  lateral  acceleration  is 
split  equally  between  the  aircraft  one  can  find  the  new  value 
for  yj.  With  these  values  the  conflict  geometry  is  defined. 

The  conflict  geometry  is  valid  under  the  following  conditions: 

1)  x  is  between  zero  and  xq  to  insure  that  the  aircraft 
will  pass  in  the  future, 

2)  The  absolute  value  of  y^  is  less  than  200  because 
this  is  the  highest  observed  lateral  velocity,  and 

3)  The  absolute  value  of  y  is  less  than  400  because  this 
would  be  the  maximum  relative  lateral  velocity. 

The  procedure  outlined  above  is  applied  to  a  set  of  initial 
separations  and  closing  speeds.  These  separations  and  closing 
speeds  are  chosen  in  the  following  manner.  The  Initial 
alongtrack  separation  is  the  greatest  alongtrack  separation  at 
which  there  is  still  a  possibility  of  a  conflict.  This 
alongtrack  separation  is  given  as 

xQ  -  (Sx  +  S2)  *  T+D  (D-9) 

e 

where  and  S2  are  the  forward  speeds  of  the  two  aircraft, 

T  is  the  look  ahead  time  and  D  is  the  radar  separation  minimum. 
The  crosstrack  separation  ranges  from  -25  to  -1-25  nmi.  The  minus 
values  indicate  that  the  aircraft  have  "switched  positions”  with 
respect  ot  their  assigned  routes.  The  crosstrack  closing  speeds 
range  from  +400  kts  to  -400  kts.  The  minus  closing  speed 
indicates  an  opening  situation  for  positive  values  of  crosstrack 
separation  and  a  closing  situation  for  negative  values  of 
crosstrack  separation.  Therefore,  given  an  aircraft  pair  which 
is  separated  alongtrack  by  x0,  crosstrack  by  yQ,  having  a 
crosstrack  closing  speed  of  y0  while  one  aircraft  has  a 
crosstrack  speed  of  y^,  and  having  a  relative  lateral 
acceleration,  one  can  determine  where  on  the  conflict  surface, 
if  at  all,  the  trajectory  will  penetrate. 
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D.3  Detection  Delays 
P«3.1  Methodology 


The  objective  of  the  detection  delay  analysis  is  to  estimate  a 
distribution  or  distributions  of  the  detection  delay.  This 
delay  is  the  time  it  takes  to  recognize  the  penetration  of  a 
trajectory  which  is  really  inside  the  conflict  surface*  This 
estimate  was  made  via  simulation.  The  flow  of  the  simulation 
proceeded  as  follows: 

1)  Randomly  choose  an  initial  geometry  (i.e.  x**35  nmi,  y, 
y,  y^  and  y)  and  find  where  the  trajectory  penetrates  the 
conflict  surface  as  previously  described, 

2)  Select  a  surveillance/tracker  error  from  a  correlated 
quadri variate  normal  distribution  (see  Appendix  B), 

3)  Add  the  position  and  velocity  errors  due  to  the 
surveillance/tracker  system  to  the  nominal  positions  and 
velocities  and  check  to  see  whether  the  observed  positions 
and  velocities  place  the  aircraft  pair  outside  the  conflict 
surface, 

4)  Repeat  steps  2)  and  3)  100  times  and  estimate  the 
probability  of  not  detecting  the  pair  within  the  conflict 
surface , 

5)  Estimate  the  probability  of  first  detecting  the  pair  of 
aircraft  within  the  conflict  surface  by  taking  the  product 
of  the  probability  of  detecting  the  pair  inside  the 
conflict  surface  on  this  update  and  the  probability  of  not 
detecting  the  pair  inside  the  conflict  surface  on  all  the 
previous  updates, 

6)  Project  the  trajectory  of  the  aircraft  pair  ahead  over 
the  update  time  interval, 

7)  Repeat  steps  2)  through  6)  either  10  times  or  until  the 
probability  of  first  detection  at  the  current  update  is 
less  than  10"7, 


8)  Repeat  steps  1)  through  7)  for  the  given  number  of 
replications. 


The  generation  of  the  survelllance/track  errors  is  based  on 
sampling  from  a  correlated  quadrivariate  normal  distribution. 

The  sampling  procedure  is  explained  in  Appendix  B.  The  four 
variables  of  this  distribution  are  the  alongtrack  separation,  x, 
the  crosstrack  separation,  y,  the  alongtrack  closing  speed,  x, 
and  the  crosstrack  closing  speed,  y.  The  required 
transformations  to  generate  the  required  variates  are 

x  -  (0.28151^  +  0.0821U2  -  0.02U3  -  0.06U4)  <rx 

x  -  (53.6073 U2  +  5.5U3  -  11.0U4)  or* 

y  -  (0.12U3  +  0.16Uy)  <ry 

y  -  (36.0U4)  o" y 

where  U3,  U2,  U3,  and  U4  are  four  independent  normally 
distributed  [N(0,1)J  variables.  The  values  of  the  sigmas  are 
given  in  Appendix  B. 

The  determination  of  whether  a  particular  geometry  (l.e.  set  of 
x,  y.  y,  i)  is  inside  or  outside  the  conflict  surface  is  made  by 
basically  using  the  conflict  alert  algorithm. 

D.3.2  Results 

The  simulation  was  run  for  10  penetration  points.  The 
penetration  points  were  randomly  chosen  with  two  basic 
constraints.  These  constraints  were  that  some  penetrations  must 
occur  on  the  leading  edge  of  the  conflict  surface  and  some  of 
the  penetrations  must  have  significant  accelerations.  The  rest 
of  the  penetrations  would  be  backside  penetrations  and/or  have 
small  accelerations  (  <  >015g  combined  lateral  acceleration 
between  the  two  aircraft).  These  constraints  are  based  on  two 
hypotheses:  1;  that  it  may  be  the  backside  or  leading  edge 
penetration  position  that  determines  the  detection  delay  and  2) 
that  it  may  be  the  magnitude  of  the  acceleration  that  determines 
the  detection  delay. 

The  results  on  Table  D-l  show  that  there  are  4  leading  edge 
penetrations  (the  "LN  in  the  LBI  column)  ar.  I  6  backside 
penetrations.  The  cumulative  hi'tograms  of  the  probability  of 
first  detection  for  each  penetration  are  shown  in  Table  D-2. 
These  cumulative  histograms  are  plotted  in  Figure  D-l.  As  can 
be  seen  the  leading  edge  distribution  tends  to  rise  faster  than 
the  backside  distribution.  The  reason  for  this  is  conjectured 
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FIGURE  D-1 

CUMULATIVE  DETECTION  DELAY  DISTRIBUTIONS 


to  be  the  following.  The  leading  edge  of  the  conflict  surface 
Is  curved  and  thus  for  an  aircraft  pair  on  the  surface  it  is 
more  likely  that  it  appear  to  be  outside  the  surface.  However, 
after  one  or  two  radar  updates  it  is  definitely  inside  the 
surface  and  the  chances  of  being  detected  there  quickly  become 
very  high.  On  the  other  hand,  an  aircraft  pair  on  the  backside 
with  not  much  acceleration  would  "slide"  down  the  backside  of 
the  conflict  surface.  At  each  update  there  would  be 
approximately  a  50-50  chance  of  detection.  Even  with 
acceleration  it  will  take  several  updates  for  the  aircraft  pair 
to  get  far  enough  into  the  backside  of  the  conflict  volume  to 
raise  the  probability  of  detection  above  0.5  on  a  given  update. 

Based  on  the  reasoning  above  it  was  decided  to  assign  two 
detection  delays  distributions.  One  is  assigned  for  aircraft 
pairs  entering  the  leading  edge  of  the  conflict  surface  and  one 
for  aircraft  pairs  entering  the  backside  of  the  conflict 
surface.  The  determination  of  the  detection  delay  distribution 
was  made  by  taking  the  average  cumulative  probability  at  each 
update  for  the  leading  edge  results  and  for  the  backside 
results.  The  resulting  histograms  out  to  120  seconds  of  delay 
are  shown  in  Table  D-3  of  the  Appendix  and  in  Figure  2-8  in 
Volume  I  of  this  report. 
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APPENDIX  E 


CONFLICT  RESOLUTION 


E.l  Resolution  Algorithm 

In  the  previously  modelled  resolution  maneuver  for  the  same 
direction  case,  one  of  the  aircraft  was  assumed  to  turn  back 
toward  its  assigned  route.  The  maximum  turn  was  such  that  both 
aircraft  were  heading  in  the  same  direction.  It  was  assumed 
that  after  such  a  turn,  the  conflict  would  be  over. 

In  the  opposite  direction  case,  it  is  not  necessarily  the  best 
maneuver  to  turn  one  of  the  aircraft  back  to  its  assigned 
route.  The  closing  speeds  are  fast  enough  that  in  some  cases  it 
is  advantageous  to  turn  away  from  one's  route  in  order  to  effect 
a  better  or  even  a  feasible  resolution.  Once  the  direction  of 
turn  is  decided,  there  is  still  the  issue  of  how  much  to  turn 
the  aircraft  to  resolve  the  conflict.  Surely  it  is  not 
necessary  for  the  two  aircraft  to  fly  parallel  courses  to 
resolve  the  conflict  when  they  originally  had  opposite  courses. 

The  way  the  opposite  direction  resolution  maneuver  will  be 
handled  will  be  to  use  the  experimental  en  route  Conflict 
Resolution  Advisory  algorithm  being  developed  for  the  FAA.  A 
detailed  description  of  the  full  algorithm  can  be  found  in 
Reference  4. 

Based  on  the  positions  and  the  velocities,  the  aircraft  are 
projected  ahead  a  given  number  of  seconds  (42  seconds  is  the 
current  parameter)  to  account  for  delays  and  the  instantaneous 
turn  to  be  made  by  the  maneuvering  aircraft.  After  the  delay,  a 
wedge  of  prohibited  relative  velocity  headings  is  computed  for 
the  faster  aircraft.  This  wedge  is  defined  by  the  vectors  UL 
and  UR  as  shown  in  Figure  E-l.  If  the  relative  velocity  vector 
is  inside  this  wedge  and  the  turn  is  made  instantaneously,  to 
take  the  relative  velocity  vector  outside  the  wedge,  the 
aircraft  will  miss  by  a  distance  SPR  (-  5  nmi). 

Now  consider  that  the  faster  aircraft  will  be  maneuvered.  The 
slower  aircraft  will  have  an  observed  track  heading.  However, 
depending  on  the  type  of  tracking  (Free  or  Flat),  there  will  be 
errors  in  the  heading  of  the  slower,  nonturning  aircraft.  This 
is  described  by  a  wedge  of  headings  V2L  to  V2R  as  shown  in 
Figure  E-2  (Aircraft  2  is  the  nonturning  aircraft  by 
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convention.)  The  problem  is  to  find  the  headings  for  the  faster 
aircraft  such  that  the  relative  velocity  vectors  do  not  fall  in 
the  wedge  defined  by  the  vectors  UL  and  UR.  Figure  E-2  depicts 
the  geometry.  If  the  heading  of  the  faster  aircraft  is  to  the 
left  of  VlL  or  to  the  right  of  V1R  no  maneuver  is  required  to 
miss  the  slower  aircraft  by  SPR.  In  the  algorithm,  it  is  also 
assumed  that  the  observed  heading  of  the  turning  aircraft  may 
also  be  in  error  so  that  a  given  number  of  degrees  is  added  to 
VlL  and  V2R  to  account  for  this  error.  If  the  heading  of  the 
faster  aircraft  were  between  VlL  and  V1R  then  it  would  have  to 
turn  to  the  left  or  right  to  get  its  velocity  vector  out  of  this 
wedge . 

If  one  were  to  turn  the  slower  aircraft  the  velocity  geometry 
would  be  as  shown  in  Figure  E-3.  In  this  particular  case,  UL, 
UR,  V2L ,  and  V2R  are  such  that  there  two  prohibited  wedges  Tor 
aircraft  1,  the  slower,  turning  aircraft.  There  could  be  no 
prohibited  wedge,  one  wedge,  two  wedges,  or  no  permitted 
velocity  vector  (i.e.,  the  aircraft  will  not  be  separated  by  SPR 
no  matter  what  the  slower  aircraft  does).  Thus,  if  the  slower 
aircraft's  velocity  vector  is  inside  a  prohibited  wedge  then  it 
must  turn  to  get  its  velocity  vector  out  of  the  wedge.  Here 
again  the  prohibited  wedges  are  augmented  to  account  for  the 
heading  errors  of  the  turning  aircraft. 

The  logic  of  the  Conflict  Monitoring  Analysis  will  choose  the 
following  maneuver.  If  neither  aircraft's  velocity  vector  is  in 
a  prohibtied  wedge,  then  neither  is  required  to  make  a  turn. 
Otherwise  the  aircraft  which  has  to  make  the  smallest  turn  will 
be  the  aircraft  designated  to  make  the  maneuver.  It  is  usually 
the  faster  aircraft  which  will  make  the  maneuver.  Thus,  both 
the  magnitude  and  direction  of  the  turn  is  determined  by  this 
algorithm. 

If  the  aircraft  are  nearly  the  same  speed,  one  cannot  be  sure 
which  aircraft  is  the  faster  one.  In  this  case,  each  aircraft 
is  forced,  in  turn,  to  be  a  certain  percentage  faster  than  the 
other.  The  slower  aircraft  is  maneuvered  in  each  case  because 
this  is  the  more  conservative  approach  (i.e.,  the  slower 
aircraft  will  have  to  turn  through  a  larger  angle).  The 
aircraft  which  is  required  to  turn  through  the  smallest  angle 
when  it  is  forced  to  be  the  slower  aircraft  is  chosen  to  be  the 
turning  aircraft. 


Prohibited  Heading 
Wedges  Aircraft  1 


FIGURE  E-3 

PROHIBITED  HEADING  WEDGES 
SLOWER  AIRCRAFT 


The  resolution  advisory  algorithm  accounts  for  the  uncertainties 
in  the  heading  errors  of  the  aircraft  due  to  the 
surveillance/tracker  errors  by  placing  wedges  about  the  measured 
headings  of  the  aircraft*  The  output  from  the  resolution 
advisory  algorithm  tells  the  controller  how  far  left  and  how  far 
right  one  of  the  aircraft  must  turn  to  get  out  of  the  conflict 
situation.  In  the  Conflict  Monitoring  Analysis  it  is  assumed 
that  the  controller  will  pick  the  direction  of  turn  to  minimize 
the  angle  through  which  the  aircraft  must  turn. 

Assume  that  the  true  initial  configuration  of  the  conflict 
geometry  is  known.  If  this  information  is  used  in  the  Conflict 
Resolution  Advisory  algorithm,  a  decision  to  turn  one  of  the 
aircraft  right  or  left  is  returned.  Now,  if  it  is  assumed  that 
the  surveillance/tracker  system  estimated  the  positions  and 
velocities  of  the  aircraft  to  have  values  different  than  the 
true  configuration  case,  it  is  possible  that  the  decision  from 
the  conflict  resolution  algorithm  could  be  to  turn  the  aircraft 
in  the  opposite  sense.  It  is  of  Interest  to  the  Conflict 
Monitoring  Analysis  to  be  able  to  estimate  the  probability  of  an 
Improper  turn  sense  by  knowing  something  about  the  true  geometry 
of  the  conflict. 

The  object  of  this  investigation  is  to  develop  an  estimate  of 
the  probability  of  advising  the  Improper  turn  sense  as  a 
function  of  the  conflict  geometry. 

Characterizing  the  conflict  geometry  is  best  done  by  considering 
the  logic  of  the  Conflict  Resolution  Advisory  algorithm.  That 
logic  is  based  on  a  wedge  of  prohibited  headings  and  the  heading 
of  the  aircraft  to  be  turned  (s*e  Figure  E-4).  If  the  heading 
is  nearer  to  the  right  (left)  boundary  of  prohibited  wedge,  the 
preferred  command  is  to  turn  right  (left).  It  seems  reasonable 
that  if  the  heading  were  close  to  the  wedge  boundary  there  would 
be  a  high  probability  of  consistently  choosing  the  preferred 
turn  sense.  It  also  seems  reasonable  that  if  the  heading  were 
in  the  middle  of  the  prohibited  wedge  there  would  be  a  50-50 
chance  that  the  advisory  would  recommend  either  a  left  or  right 
turn.  Based  on  this  argument,  a  variable  is  defined  which  is 
the  ratio  of  the  minimum  angle  between  the  aircraft  heading  and 
the  boundary  of  the  prohibited  wedge  to  the  total  angle  of  the 
prohibited  wedge.  The  value  of  this  ratio  will  vary  from  0  to 
0.5. 
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The  second  pare  of  this  problem  Is  to  estimate  Che  probability 
of  selecting  che  preferred  turn  sense  as  a  function  of  the  above 
defined  ratio.  This  estimate  is  based  on  a  simulation.  The 
simulation  was  structured  as  follows.  Initial  geometries  were 
chosen  at  random.  These  geometries  approximated  conflict 
geometries  in  that  it  was  guaranteed  that  the  time  to  closest 
approach  was  less  than  2  minutes  and  the  minimum  separation  was 
less  than  5  nml.  These  geometries  represented  the  “true” 
geometries  of  the  aircraft  pairs.  Twenty  such  geometries  were 
chosen  and  the  Conflict  Resolution  Advisory  was  applied  to 
determine  the  direction  of  the  turn  and  the  value  of  the  ratio. 
For  each  "true'  geometry  a  set  of  position  and  velocity  errors 
were  added  to  account  of  the  surveillance  and  tracking  system 
(see  Appendix  B).  These  new  position  and  velocity  estimates  can 
be  regarded  as  the  "perceived"  geometries.  For  each  of  the 
"perceived"  geometries  (there  were  50  for  each  "true"  geometry), 
the  Conflict  Resolution  Advisory  was  used  to  determine  the 
direction  of  the  turn.  The  probability  that  the  turn  based  on 
the  "perceived"  geometries  were  in  the  same  direction  as  the 
turns  based  on  the  "true"  geometry  was  estimated  as  the 
percentage  of  turns  in  the  same  direction  as  that  advised  for 
the  "true"  geometry. 

It  should  be  stated  here  that  there  are  four  basic  outcomes  from 
the  Conflict  Resolution  Advisory  When  dealing  with  the 
"perceived"  geometries.  They  are: 

1.  The  same  aircraft  gets  the  same  direction  turn. 

2.  The  same  aircraft  gets  the  opposite  direction  turn. 

3.  The  other  aircraft  gets  a  turn. 

4.  No  aircraft  gets  a  turn. 

It  was  determined  from  the  simulation  that  in  almost  all  cases, 
the  faster  aircraft  will  get  the  turn  advisory.  It  is  also 
extremely  rare  that  the  other  aircraft  would  get  an  advisory 
even  due  to  surveillance /tracker  errors.  Thus  outcome  3  is  of 
no  concern.  Outcome  4  is  also  not  of  concern  by  definition.  It 
is  assumed  that  the  aircraft  pair  is  within  the  conflict 
boundary  so  that  this  necessitates  a  turn  based  on  the  Conflict 
Resolution  Advisory  logic.  If  the  "perceived”  geometry  leads  to 
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a  no  turn  outcome,  then  It  necessarily  follows  that  the  aircraft 
pair  has  not  yet  been  perceived  within  the  conflict  boundary. 

The  delay  in  this  perception  is  accounted  for  in  the  detection 
delay  analysis  discussed  in  Appendix  D. 

Thus  the  conditional  probability  that  the  same  aircraft  receive 
either  the  same  or  opposite  direction  turns  given  that  a  turn  is 
recommended  is  of  interest.  The  probability  of  a  proper 
direction  turn  is  shown  in  Figure  E-5.  This  probability  is 
plotted  against  the  ratio  defined  above.  The  points  represent 
outcomes  of  the  simulation.  The  line  is  the  regression  fit  to 
the  points.  The  regression  line  is 

y  -  1.0888  Jo.25  -  x2  +  0.45009 

where  y  is  probability  of  the  proper  turn  sense  and  x  is  the 
ratio  of  the  minimum  angle  between  the  aircraft  heading  and  the 
boundary  of  the  prohibited  wedge  to  the  total  angle  of  the 
prohibited  wedge.  The  R^  value  of  this  regression  is  0.944. 

As  one  can  see  from  Figure  E-5  with  the  heading  near  to  the  edge 
of  the  prohibited  wedge  (Ratio  near  to  zero)  the  probability  is 
high  that  the  surveillance/t.acker  error  will  not  cause  an 
Improper  turn  sense  to  be  recommended.  However,  if  the  heading 
is  near  the  middle  of  the  prohibited  wedge  the 
surveillance/tracker  errors  are  as  likely  to  advise  a  turn  in 
one  direction  as  In  the  other. 


APPENDIX  F 


OPPOSITE  DIRECTION  OVERLAP  REGIONS 


This  appendix  documents  the  equations  and  methods  used  to  compute 
the  overlap  regions  for  opposite  direction  traffic.  The  basic 
difference  between  this  work  and  that  done  for  the  same  direction 
traffic  is  that  accelerations  are  considered.  This,  in  turn, 
introduces  non-linearities  into  the  equations  that  were  not  in  the 
same  direction  case. 

The  problem,  simply  stated,  is  to  find  the  combinations  of  delay 
time  and  turn  rate  of  one  of  the  aircraft  such  that  an  aircraft  pair 
starting  with  a  particular  geometry  will  end  in  an  overlap 
condition  during  the  course  of  a  single  maneuver.  If  the  aircraft 
are  represented  as  two  circles  each  of  radius  R,  then 

(x2(t)  -  x1(t))2+(y2(t)-y1(t))2  <(2R)2  (F-l) 

is  the  condition  of  overlap  after  t  seconds.  The  length  of  time  t 
is  composed  of  delay  time,  t^,  and  the  time  in  the  turn,  tt. 

Thus 


c  "  cd  +  ct 


(F-2) 


The  time  in  the  turn  depends  on  the  angle  through  which  the  aircraft 
turns  and  the  turn  rate.  The  angle  is  denoted  as  P 1 17 1  where  T]  is 
the  angle  from  the  Conflict  Resolution  Advisory  algorithm  (Appendix 
E)  and  P  is  the  fraction  of  the  turn  completed  before  overlap.  P 
ranges  from  0  to  1.  The  sign  convention  on  T]  is  positive  for  left 
turns  and  negative  for  right  turns.  The  turn  rate  is  denoted  as  <J. 
The  sign  convention  for  u)  is  the  same  as  that  of  77.  Thus: 


There  are  certain  assumptions  used  in  specifying  the  acceleration. 
First  the  acceleration  which  is  specified  is  the  relative  average 
crosstrack  acceleration  between  the  two  aircraft.  Furthermore,  this 
crosstrack  acceleration  is  assumed  to  be  constant  and  divided 
equally  between  the  two  aircraft.  Another  basic  assumption  in  the 
Conflict  Monitoring  Analysis  is  that  the  forward  speeds  of  the 
aircraft  are  also  constant.  It  should  be  recognized  that  a  constant 
lateral  acceleration  and  a  constant  forward  speed  are  inconsistent 
assumptions.  What  was  done  was  to  adjust  the  alongtrack  track 
velocity  to  make  the  forward  velocity  nearly  constant  over  the 


F-l 


length  of  Che  encounter.  It  .should  be  noted  that  the  resulting 
trajectory  will  not  be  circular  because  of  the  constant  crosstrack 
acceleration. 


The  expression  that  Is  used  for  the  alongtrack  velocity,  for 
aircraft  :  Is 


x,  ■  x  . 
1  ol 


-  V  <*ol  +  V  > 


(F-4) 


ol 


where 


xoi 

*ol 

*i 


the  Initial  alongtrack  velocity 
the  initial  crosstrack  velocity 


-  the  crosstrack  acceleration  for  that  aircraft 
(one-half  the  relative  crosstrack  acceleration). 


When  an  aircraft  turns  its  alongtrack  velocity  will  be  changed  by 


—  (cos£  -  cos  (  £+Wtt)) 


and  its  crosstrack  velocity  will  be  changed  by 
(sin  (  £  +  wt  )  -  sin  $  ) 

U)  t 


(F-5; 

(F-6; 


where 


V  ■  the  forward  speed  of  the  turning  aircraft 


£  -  the  angle  between  the  instantaneous  velocity  vector  at 
the  start  of  the  turn  and  the  positive  crosstrack 
axis,  (see  Figure  F-l) 


The  definition  of  £  is  complicated  by  the  acceleration  and  the  fact 
that  it  should  to  be  in  terms  of  the  initial  parameters.  After  a 
time  delay  tg,  the  angle  4  *8  given  by 


kfi 


r  v-  Vd  V  Vdj 

«(t")-“0"1 - r-rr? - 

yoi  yi  cd 


(F-7) 


Conbining  (F-2),  (F-3),  (F-4),  (F-5)  and  (F-6)  we  arrive  at 


x„-  x, »  x  t.i  -  2  A  t.t  -  Bt  .t?“(A~D)t j  -  Btj  t 
2  1  a  d  t  at  a  at 

-(B-C)t^  +  Ett-  At^  -  B  /3 
+  Vtt  f"  cos  £  -  cos  (  £  +  P|f?l)l 

D  I  I  L 


(F-8) 


where 


y2‘yi“  y+tdy  +y02tt+(y2-y1)  +  y2tdtt 

2 

+  [81“  U 


A-  y2yo2 


ol 

D-  *l*ol 


(F-9) 


E-  *o2 

y0i  -  Initial  crosstrack  speed  of  aircraft  1 
xoi  *  Initial  alongtrack  speed  of  aircraft  1 
Yl  •  crosstrack  acceleration  of  aircraft  1 


x  -  Xo2  -  Xoi 

y  ■  y02  -  yQi 


The  top  sign  is  appropriate  when  U)  >  0  while  the  bottom  sign  is 
appropriate  when  u  <  0.  It  should  be  noted  that  the  definitions 
of  the  variables  are  based  on  the  cartesian  coordinate  system 
shown  in  Figure  F-l.  This  convention  varies  from  that  used  to 
define  the  variables  in  Appendix  C. 

The  ultimate  aim  is  to  find  those  values  of  td  and  (j  which 
satisfy  (F-l).  This  can  be  done  by  considering  that  (F-l)  has 
the  general  form  of  the  circle.  Therefore  (F-l)  can  be 
partitioned  into  two  equations  by  defining  a  new  variable. 


x2  “  xi  -  2R  cos  6  -  0 
y2  -  yi  -  2R  sin  6=0 


0  <  0  $  2tt 


(F-10) 


There  are  now  two  equations  and  two  unknowns.  However,  (F-10) 
are  nonlinear  equations.  The  solution  to  these  equations  is 
found  by  using  a  Newton's  iteration  technique  to  solve  for  td 
and  tt.  The  value  of  G)  is  then  found  as 


As  an  example,  six  values  of  0  can  be  chosen.  This  will  yield  a 
region  of  overlap  in  td  -  <j  space.  One  such  overlap  region 
is  shown  in  Figure  F-2.  If  one  computes  the  distance  between 
the  aircraft  as  a  function  of  time  for  the  case  of  td  ■  13.88 
sec  followed  by  a  turn  at  .19575  deg/sec  one  finds  that  the 
aircraft  come  within  a  distance  2R  at  about  125.7  seconds  (see 
Figure  F-3). 

Given  ranges  of  x,  y,  y,  yj.,  and  y,  many  of  these  overlap 
regions  can  be  computed.  The  result  of  doing  this  for  one 
particular  set  of  such  ranges  of  values  is  shown  in  Figure  2-12 
in  Volume  I  of  this  report. 


<0.19575,  13.88) 


x_  “  31.5  nmi 
yi  “  -1.67  oni 
i,  -  -389.5  kt8 
yt  -  -20.0  kts 
y2  -  0  kts/hr 


1 - 

0.195 


R  “  0.0126  nml 


0.196 


1 - 

0.197 


TURN  RATE  (DEG/SEC) 


FIGURE  F-2 

EXAMPLE  OF  OVERLAP  REGION 


SEPARATION  (NMI) 


RQURE  F-3 

SEPARATION  AS  A  FUNCTION  OF  TIME 
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APPENDIX  G 


PROBABILITY  OF  HORIZONTAL  OVERLAP 

In  Section  2  of  Volume  I  of  this  report  the  procedure  for  estimating 
the  probability  of  horizontal  overlap  was  outlined*  There  was  one 
term  In  that  computational  process  which  is  the  probability  of 
getting  into  overlap  given  that  the  aircraft  pair  started  in 
conflict  surface  cell  j*  The  estimation  of  this  probability  is 
discussed  in  this  appendix. 

The  estimate  of  the  probability  of  overlap  given  the  aircraft  pair 
started  in  cell  j  is  based  on  the  following  expression: 

P(overlap  |  CSj)  -  (P£ *Pp+(l-P£ )*(l-Pp))*P£tdo  (6-1) 

+((l-Pf)*Pp+P£(l-Pp»*Prtdo 

where 

Pj,  *  probability  of  left  turn  command 

Pp  -  probability  of  proper  sense  turn 

pftdo  *  probability  of  left  turn,  delay,  and  overlap 

region 

prtdo  ■  probability  of  right  turn,  delay,  and 

overlap  region 

Each  of  these  probabilities  will  now  be  discussed. 

From  Volume  I  of  this  report  recall  that  from  a  cell  on  the  conflict 
surface  an  overlap  region  in  delay/turn  rate  space  can  be  created. 

An  example  of  such  an  overlap  region  is  given  in  Figure  G-l.  For 
the  particular  region  shown  both  left  (positive  turn  rate  )  and 
right  (negative  turn  rate  )  turns  could  lead  to  overlap.  In  most 
conflict  surface  cells  only  one  direction  turn  will  lead  to  an 
overlap.  To  estimate  the  probabilities  P£tdo  aiu*  prtdo  one  luis 
to  consider  the  turn  rate  probability  distribution  and  the  total 
delay  probability  distribution.  One  then  computes  the  probability 
of  having  a  range  of  turn  rates  and  a  range  of  delay  and  being  in 
the  overlap  region.  A  discussion  of  the  computation  of  this 
probability  is  found  in  Appendix  B  of  Volume  II  of  Reference  3. 

The  probability  of  a  left  turn  command  is  made  by  sampling  over  the 
conflict  surface  cell.  At  each  point  the  Conflict  Resolution 
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Advisory  algorithm  is  used  to  indicate  whether  a  right  or  left  turn 
should  be  made.  The  value  P  is  estimated  by  the  fraction  of  left 
turn  commands  from  the  sample.  In  most  cases  for  a  given  conflict 
surface  cell  the  value  of  P  will  be  nearly  1  or  nearly  0.  It  is 
also  true  that  in  most  cases  the  predominant  turn  direction  is  away 
from  the  overlap  situation.  This  indicates  that  with  good 
surveillance  information  the  proper  turn  direction  is  given  by  the 
Conflict  Resolution  Advisory. 

The  term  Pp  is  the  probability  that  the  proper  turn  direction  is 
in  fact  given  to  the  pilot.  An  improper  turn  direction  would  be  due 
to  surveillance/tracker  error.  The  estimate  of  Pp  is  based  on  the 
heading  of  the  turning  aircraft  with  respect  to  the  wedge  of 
prohibited  headings  given  in  the  Conflict  Resolution  Advisory 
algorithm  (see  Appendix  E) .  From  a  simulation  it  was  found  that  if 
the  aircraft  heading  were  near  the  middle  of  the  prohibited  wedge  of 
headings  (see  Figure  6-2)  the  chances  are  50-50  that  the  improper 
turn  sense  would  be  selected  by  the  Conflict  Resolution  Advisory 
algorithm.  If  the  aircraft's  heading  were  near  the  edge  of  the 
prohibited  wedge  of  headings  the  chances  are  remote  that  the 
surveillance  errors  would  cause  an  improper  turn  sense.  A 
simulation  was  performed  to  arrive  at  an  estimate  of  the  probability 
of  the  Conflict  Resolution  Advisory  algorithm  selecting  the  proper 
turn  sense  as  a  function  of  the  position  of  the  aircraft  heading 
within  the  wedge  of  prohibited  headings.  This  simulation  is 
discussed  in  Appendix  E.  The  result  is  that 

Pp  -  1.0888  70.25  -  R2  +  0.45009 

where  R  is  the  ratio  of  the  angle  through  which  the  aircraft  is  to 
turn  to  the  total  angle  covered  by  the  prohibited  wedge. 
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